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ABSTRACT 

We study the dipole structure of an extended redshift sample of Abell/ACO 
clusters in order to infer the dynamical origin of the motion of the Local Group 
(LG). To further elucidate the constraints this motion places on dark matter 
models, we use numerical simulations based on an optimized version of the 
truncated Zel'dovich approximation which we have shown in previous work to 
provide a reliable representation of large-scale gravitational clustering. Taking 
advantage of their low computational cost, we run 20 realizations of each of six 
different dark matter models: four of these have a density parameter fl = 1, 
while the other two have fl = 0.2, one with and one without a cosmological 
constant term. For the Abell/ACO sample, we have evaluated the parameter 
P = f2°' 6 /6 c ; (where b c \ is the linear bias parameter for the clusters), which 
reaches its asymptotic value at R CO nv — 160 h~ x Mpc. Convergence occurs at 
this scale whether calculations are performed in the LG or CMB frame, but 
the asymptotic value differs in these two cases: we find Plg = 0.15 ± 0.04 and 
Pcmb = 0.25 ± 0.06, respectively. After identifying in the simulations those 
observers having local densities and peculiar velocities similar to those of the 
Local Group, we construct mock cluster samples around them reproducing the 
same observational biases, and apply to these mock samples the same method 
of analysis as we used for the Abell/ACO sample. We find that an alignment 
between the cluster dipole and observer velocity ('CMB' dipole) directions, such 
as that observed (A6^ 20°), should not be expected necessarily: much larger 
misalignment angles are often found in all models considered. This, together 
with the large observer-to-observer variance estimates of /3, makes it difficult to 
place any firm constraints on cosmological models. This result suggests that the 
dipole analysis of the cluster distribution has a relevance that is cosmographical, 
rather than cosmological. Our results also demonstrate that the large amplitude 
and convergence depth of the observed cluster dipole cannot be taken as strong 
evidence either for or against a low-density Universe. 

1 INTRODUCTION provide a means to determine the cosmological density 

parameter Q . According to the linear theory of gravi- 
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ity of an observer placed at the centre is given in terms 
of the density fluctuation field S by 



H f(U a ) f 



where 



S(r) 



P( r ) 



(1) 



(2) 



In eq.(l), /(f2 ) — 6 to a good accuracy, virtu- 
ally independently of the value of the cosmological con- 
stant (e.g. Peebles 1993); the window function W filters 
out small-scale fluctuations and defines the "size" and 
"shape" of the observer through its functional depen- 
dence on r. As the size of the sphere is increased, the 
contribution to v v increases until it eventually saturates 
at scales where f < 1. 

Since the only well-defined cosmic peculiar velocity 
is that of the Local Group of galaxies (LG hereafter), 
determined from the dipole pattern of the CMB radi- 
ation temperature (Kogut et al. 1993), the relation (1) 
is most naturally applied to the LG. By measuring the 
integral in eq.(l), we could obtain an estimate of Q . 
Although easy in principle, it is a hard task in practice, 
since we have no direct way to measure 6(r). The best 
we can hope to do is to use some class of cosmic objects 
as tracers of the underlying density field and hope that 
their point -like distribution, n(r) = Y2i^ D ( r ~ r ') * s 
simply related to the underlying 6(r). In the previous 
expression Sb stands for the Dirac delta-function. The 
assumption usually invoked is that relative fluctuations 
in the object number counts and matter density fluc- 
tuations are proportional to each other, at least within 
sufficiently large volumes, according to the linear bias- 
ing prescription (cf. Kaiser 1984; Bardeen et al. 1986): 



6n(r) 



(3) 



n p 

where b is what is usually called biasing parameter. If 
this is the case, then the comparison between velocity 
and acceleration provides a constraint on the parameter 
(3, given by 

e.lfl^. (4 ) 

It should be clear that in order to infer the /3-parameter 
by relating the dipole estimated from an observational 
sample to the LG velocity, one must be sure that (a) the 
sample size is sufficiently large to contain all the fluctu- 
ations responsible for the LG motion, and (fe) the CMB 
and the mass tracer dipoles are sufficiently well aligned 
to ensure that linear theory, as expressed by eq.(l), is 
applicable and that their estimates are not contami- 
nated by the noise associated to the sparse sampling 
of the continuous density field. 



The dipoles of various populations of extragalac- 
tic objects have been determined up to this date: opti- 
cal galaxies (Lahav 1987; Plionis 1988; Lahav, Rowan- 
Robinson & Lynden-Bell 1988; Lynden-Bell, Lahav & 
Burstein 1989; Hudson 1993), IRAS galaxies (Meiksin 
& Davis 1986; Yahil, Walker & Rowan-Robinson 1986; 
Villumsen & Strauss 1987; Strauss & Davis 1988; 
Rowan-Robinson et al. 1990; Strauss et al. 1992; Plio- 
nis, Coles & Catelan 1993); X-ray active galactic nu- 
clei (Miyaji & Boldt 1990); X-ray clusters (Lahav et 
al. 1989) and Abell clusters (Scaramella, Vettolani & 
Zamorani 1991; Plionis & Valdarnini 1991 [hereafter 
PV91]; Branchini & Plionis 1995). In all cases the dipole 
moment is reasonably well aligned with the CMB dipole 
suggesting that gravity is indeed responsible for the Lo- 
cal Group motion and that linear theory (eq. 1) can be 
applied to infer the value of (3. However, these stud- 
ies have provided conflicting estimates of (3 which, in 
some cases, could be attributed to the different amount 
of biasing with respect to the background matter dis- 
tribution (i.e. different b values) but, in other cases, 
one would have to infer differences also in the value 
of Q (for an overview, see Dekel 1994). In particular, 
IRAS galaxies have provided values /3iras ~ 0.8 ± 0.2 
implying (0 ,biRAs) ~ (1,1.2), optical galaxies val- 
ues somewhat smaller of f3 op t ~ 0.6 ± 0.2, implying 
(Q ,b op t) ~ (1,1.7) and clusters even smaller values 
f3 c i ~ 0.15 ± 0.05 implying either a quite large cluster 
biasing parameter (b c i ~ 7) which is not supported by 
correlation function and power-spectrum studies (Pea- 
cock & Dodds 1994; Jing & Valdarnini 1993) or a value 
of n < 1 (Plionis, Coles & Catelan 1993; Coles & Ellis 
1994; see Branchini & Plionis 1995 for /3 ci ~ 0.2 which 
is consistent with Q PS 1). 

In this work we plan to exploit the value of (3 de- 
rived from the cluster distribution to address the follow- 
ing two questions: 

(i) Is the use of linear theory to infer the value of f3 c i 
(and, possibly, of f2 ) from an observational setup 
similar to that of the real Abell/ ACO cluster sam- 
ple reliable? 

(ii) Does the cluster dipole provide a strong con- 
straint to discriminate between different Dark Mat- 
ter (DM) models? 

The paper is organized as follows. In Section 2 we 
explain our method of analysis. In Section 3 we describe 
the Abell/ ACO cluster sample used and the results ob- 
tained in the analysis of the real cluster dipole. In Sec- 
tion 4 we present the simulation method applied to gen- 
erate realistic cluster distributions and the 'observer' 
selection procedure. In Section 5 we discuss the results 
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of our analysis of simulated cluster catalogues. Finally, 
in Section 6, we present our main conclusions. In the 
Appendix, some technical details on the spherical har- 
monics mask method are discussed. 

2 THE DIPOLE ANALYSIS 

Under the assumption of linear gravitational instability 
and linear biasing, eq.(l) relates the observer peculiar 
velocity to the distribution of clusters according to 

N(R) 

v > = I>*h[ = p v « ■■ ( 5 ) 

2=1 

f3 = f2°' 6 /b c i; b c i is the cluster biasing parameter; N(R) 
is the total number of clusters within a distance R, Wi 
are suitable weights which account for the different clus- 
ter masses and for selection biases in the observed clus- 
ter distribution (e.g. Galactic absorption, redshift selec- 
tion, systematic biases, etc.). Estimating V c i from the 
observed cluster distribution and knowing v v from the 
CMB dipole amplitude, eq.(5) allows one to determine 
the parameter (3. 

It is useful to introduce the monopole and dipole of 
the cluster distribution evaluated at a scale R, as 

N(R) N(R) 

2=1 2=1 

According to the above expressions, if all the mass (clus- 
ters) were concentrated at a single point, then \D\ = 
ZM; for a completely isotropic distribution D = 0, while 
for a Poissonian distribution \D\ ~ ZM/\/~N. In gen- 
eral the dipole amplitude is expected to increase with 
distance for as long as the sphere of radius R has not 
encompassed the largest inhomogeneity in the cluster 
distribution. Its flattening is the signature that isotropy 
has been reached. The monopole grows as M(R) cx R 
for a uniform distribution and, therefore, if the conver- 
gence scale of the dipole R CO nv is within the sample size, 
then for scales R > R CO nv we have M(R) = 4irnR. In 
what follows we define as the effective average cluster 
density the value derived from the monopole at the scale 
of the dipole convergence: n = M{R conv )/{4TrR con - v ). 
By inserting eqs.(6) into eq.(l), we get 

which will be used in the following as the estimator 
of V c i. Note, however, that eqs.(6) can be directly ap- 
plied only in the case in which the cluster sample cov- 
ers the whole sky. This ideal situation is not satisfied 
in all the realistic cases: Galactic obscuration progres- 
sively degrades cluster detections at low Galactic lati- 



tudes and so affects the determination of both ampli- 
tude and direction of D . Different methods for correct- 
ing for this effect have been used by different authors. 
One approach is to fill the obscured portion of the sky 
uniformly (cf. Strauss & Davis 1988; Lahav, Rowan- 
Robinson & Lynden-Bell 1988). Another possibility is 
to "clone" the observed distribution near the zone of 
avoidance and extrapolate it across the Galactic plane 
(cf. Lynden-Bell, Lahav & Burstein 1989; Branchini 
& Plionis 1995). Yet an other approach is that based 
on the spherical harmonic reconstruction method. This 
method, which we will use in our present analysis, has 
been applied to the galaxy distribution by Yahil, Walker 
& Rowan-Robinson (1986), Lahav (1987), Plionis (1988, 
1989) and, for the cluster distribution, by PV91. Details 
of this method can be found in the Appendix. 

3 ANALYSIS OF THE ABELL/ACO 
CLUSTERS 

3.1 The sample 

The Abell/ACO sample (Abell 1958; Abell, Corwin & 
Olowin 1989) considered here is the same as that de- 
scribed in Borgani et al. (1995, Paper I) with additional 
new redshifts provided by the ESO Abell cluster sur- 
vey (Katgert et al. 1995). It consists of all clusters hav- 
ing both mio < 17, where mio is the magnitude of the 
tenth brightest cluster galaxy in the magnitude system 
corrected according to PV91, and an estimated distance 
< 240 ft." 1 Mpc. There are in total 261 Abell clusters 
(dec > —17°), all having measured redshifts, and 201 
ACO clusters, out of which 35 have z estimated from 
the mio-z relation derived in PV91. Since we will com- 
pare the real data with simulations based both on flat 
and open cosmological models, with and without a cos- 
mological constant, we convert redshifts into cluster dis- 
tances using the general formula for the distance by ap- 
parent size: 

r{z) = 2r o (i-n -n A ) sinh {i-n„-n A ) ' (8) 

where c is the speed of light and 

E(z) = ^a {i + z y + a A + (i - a - n A )(i + zf (9) 

(e.g. Peebles 1993). Final results are actually quite in- 
sensitive to the choice of the cosmological parameters 
to be inserted in eq.(8) so, except where differently 
specified, we will from now on use distances based on 
(f2 , fJ A ) = (0.4,0), which are intermediate between the 
cases with (fio,fi A ) = (1,0) and (fio,fi A ) = (0.2,0.8) 
cases. 
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To a good approximation, the redshift selection 
function, P(z), has the functional form (cf. Postman 
et al. 1989): 



P(z) 



f 1 

\ Ae~ z/z ° if 



Z < Z c 

z > z c 



(10) 



where A = e Zc ^ z " and z c is the redshift out to which 
the space density of clusters remains constant (volume- 
limited regime). Using a best-fitting procedure, we ob- 
tain z c = 0.0786, z = 0.01 and z c = 0.069, z a = 0.009 
for Abell and ACO samples, respectively. Since correct- 
ing for this effect (by weighting clusters by the inverse of 
P(z)) could introduce large shot noise errors, because of 
the exponential decrease of the cluster density, we pre- 
fer to limit our analysis to r ma x = 240 ft. -1 Mpc. This is 
also consistent with the box size (480 ft." 1 Mpc) of the 
simulations we shall discuss later (see below). 

In Figure 1 we present the cluster number den- 
sity estimated in equal proper volume shells and cor- 
rected for Galactic absorption, in both the Abell and 
ACO samples (with \b\ > 30°). Within the volume- 
limited regime we have riAbeii — 1.5 x 10 -5 ( ft -1 Mpc) -3 
and tiaco — 2.3 x 10 -5 ( ft -1 Mpc) -3 , respectively. 
These values correspond to average cluster separations 
of {r Abell) ~ 40 ft" 1 Mpc and (r A co) ~ 35 ft _1 Mpc. 
The difference in densities of Abell and ACO clusters 
is partly intrinsic, due to the presence of the Shap- 
ley concentration in the ACO sample (Shapley 1930; 
Scaramella et al. 1989; Raychaudhury 1989), but prob- 
ably mostly spurious, due to the higher sensitivity of 
the Illa-J emulsion plates on which the ACO sur- 
vey is based (for more details see Batuski et al. 1989; 
Scaramella et al. 1990; PV91). In order to account for 
the difference in the cluster number density we weight 
each Abell cluster according to a radial weighting func- 
tion, w(R) = uac o{R) I n Abeii{K) , which accounts for a 
distance dependence of the number density difference. 
We have also verified that our final results are insensi- 
tive in weighting the Abell or ACO clusters with w(R) 
or l/w(R), respectively. 

In order to account for cluster selection as a func- 
tion of Galactic latitude, we model the Galactic absorp- 
tion according to (see PV91): 



0(#) 



exp(-yl|cos^|- 1 ) if|6|>|6| Kr 



(11) 



if |6| < |6|iim 

with -3 = 90° — b. We calculate the parameter A by re- 
quiring eq.(ll) to fit the cluster data for \b\ > \bu m \ (= 
13°). This choice of |&;; m | is a compromise between the 
need to have a large sky coverage and, in particular, to 
contain the Perseus cluster, which was shown in PV91 
to have a significant effect on the cluster dipole, whilst 



simultaneously avoiding noise associated with regions 
of strong absorption. The best-fit values for b > 13° 
(North Galactic Cap, NGC) and b < -13° (South 
Galactic Cap, SGC) are An ^ 0.75 and As ^ 0.63, 
respectively. We use eq.(ll) for the Galactic obscura- 
tion mask to work out the 'whole-sky' monopole and 
dipole coefficients (see Appendix). 

3.2 Results 

Dipole analyses of the Abell/ ACO cluster distribution 
have been already presented by different authors (PV91; 
Scaramella, Vettolani & Zamorani 1991; Branchini & 
Plionis 1995). Here we briefly discuss our analysis of the 
data in order to make a homogeneous comparison with 
the simulation results. Note that here we do not weight 
clusters according to their mass, as done in most other 
similar studies, since we are interested in a consistent 
and meaningful comparison with results from our sim- 
ulations, in which we cannot resolve meaningfully the 
individual cluster masses. 

In Figure 2 we plot V c i evaluated according to 
eq.(7), after correcting the dipole and monopole esti- 
mates for Galactic absorption according to the spher- 
ical harmonic reconstruction method described in the 
Appendix. Upper and lower curves refer to estimates 
made in the LG and CMB frames, respectively. The 
estimate of the total dipole uncertainty includes Poisso- 
nian errors, the effects of uncertainties in the shape of 
0($) and in the radial weighting function, which PV91 
estimated to be ~ 12%. For sake of clarity, we plot the 
total uncertainty only on the dipole value at the largest 
distance. 

The overall shape of V c i(R) in both reference frames 
is quite similar: a steep rise up to ~ 40 ft" 1 Mpc reflect- 
ing the dynamical effect of the "Great Attractor" re- 
gion, with a subsequent decrease due to the pull caused 
by the Perseus-Pisces region in the opposite direction 
at ~ 50 h -1 Mpc. After that, V c i(R) remains rather flat 
up to R ~ 100 ft -1 Mpc and, then, starts increasing 
again before reaching a point of apparent convergence, 
within the completeness limit of the sample, at about 
Rconv ~ 160 ft _1 Mpc. 

Despite the similarity of the LG and CMB dipole 
shapes, their asymptotic amplitudes turn out to be 
quite different, with V" c f G = 3400 ± 500fcms _1 and 
v cmb = 2000 ± 300 fcms" 1 . This can be understood by 
noting that the distribution of clusters in redshift space 
differs from the real space (i.e. truly three-dimensional) 
distribution by the non-linear term: 



cz = H,\r\ + [v(r) - v(0)] ■ —r, 

\r\ 



(12) 
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Table 1. Abell/ACO cluster analysis. Columns 2 and 4: /3 
estimates in CMB and LG frames, after subtracting a Virgo- 
centric infall of 200 fern s~ 1 . Columns 3 and 5: misalignment 
angle between the cluster dipole and the CMB dipole in the 
two frames. Results are obtained after correcting for the 0— 
mask. 

R (h^Mpc) ficMB &-6cMB PlG A^iG 



70 


0.30 ±0.16 


24.3° 


0.19 ±0.11 


23.8° 


110 


0.29 ±0.10 


19.9° 


0.17 ±0.06 


18.5° 


150 


0.24 ±0.07 


22.3° 


0.15 ±0.05 


18.1° 


190 


0.24 ±0.06 


20.1° 


0.15 ±0.04 


17.8° 


230 


0.25 ±0.06 


22.1° 


0.15 ±0.04 


16.9° 



where v(r) is the peculiar velocity of a cluster at the po- 
sition r, in the observer frame, and v(0) is the observer 
peculiar velocity. In the LG frame (where v(0) = 620 
km/sec), structures placed in the direction of our mo- 
tion appear at a redshift smaller than that measured 
in the CMB frame (where v(0) = 0). These structures 
consequently enhance the apparent value of V c i. One is 
tempted then to prefer to make calculations in the CMB 
frame. However, if nearby clusters, which dominate the 
dipole signal, participate together with the LG in a bulk 
motion [i.e. v(r) ftj v (0)] then a better frame in which 
to evaluate the cluster dipole would be the LG frame 
since in this case we have that cz ftj H \r\. These quali- 
tative arguments serve only to indicate that one should 
consider the dipole estimated in the LG frame as an 
upper limit to the true 3-D dipole while that estimated 
in the CMB frame is a lower limit. In fact, a detailed 
3-D reconstruction procedure used to recover the real- 
space cluster distribution shows that the value of the 
3-D cluster dipole is roughly half way between the two 
above extremes (Branchini & Plionis 1995). 

Using the above values of V c i and vlg = 
620fcms _1 , we find /3 LG = 0.18 ± 0.04 and /3 C mb = 
0.33±0.06. Note, however, that these values are overesti- 
mated with respect to the real ones because a significant 
part of the LG velocity is due to the pull of the nearby 
Abell-like Virgo cluster, which is not included in the 
Abell sample, due to its proximity. Therefore, a consis- 
tent comparison between vlg and V c i can be done only 
after subtracting the contribution of Virgo from the for- 
mer. Allowing for a Virgocentric infall of ~ 180 km s -1 , 
we end up with vlg — 500fcms -1 , which brings the 
above estimates down to /3lg — 0.15 and /3cmb — 0.25. 

In Table 1 we report the values of (3 evaluated at 
different scales in the two frames, along with the mis- 
alignment angles A9 with respect to the CMB dipole 
direction. Although the final V c i convergence depth oc- 



curs at Rconv — 160 ft. -1 Mpc, its direction is already 
quite well aligned with the CMB dipole at much smaller 
distances. Furthermore, the increased coherence, due 
to redshift space distortions and/or the existence of a 
bulk flow, produces a smaller misalignment angle in the 
LG frame, with A9 LG ~ 18° while in the CMB frame 
AOcmb — 22°. 

The above values of (3 can be compared with other 
estimates of (3 based on the cluster distribution. PV91 
found from a similar analysis /3lg — 0.19 ± 0.03. Using 
their method of reconstructing real-space cluster posi- 
tions, Branchini & Plionis (1995) found (3 ~ 0.21, with 
a large uncertainty (~ 15%), while Scaramella (1995) 
found /3 LG = 0.18 ±0.02. By comparing the Abell/ACO 
cluster dipole amplitude with that of the IRAS QDOT 
galaxies, in the LG frame, Plionis (1995) estimated 
bd/bi ~ 4 for the ratio of the two biasing parame- 
ters while Peacock & Dodds (1994) suggest b c i/bi ~ 4.5 
by comparing the power-spectrum amplitude of clus- 
ters and IRAS galaxies. Together with their estimate 
of f3i = 1.0 ± 0.2 from redshift space clustering distor- 
tions, this turns into f3 c i ~ 0.22 ± 0.05 for the cluster 
distribution. All these estimates of f3 c i fall within the 
upper (CMB frame) and lower (LG frame) limits of our 
estimated /3-values. 

In the following Section, we will attempt to con- 
strain different DM models by comparing our estimate 
of (3 from the real data, with corresponding results from 
the cluster simulations. 

4 THE SIMULATIONS 
4.1 The method 

The method we use to generate realistic cluster dis- 
tributions is based on the Zel'dovich approximation 
(ZA, hereafter; Zel'dovich 1970; Shandarin & Zel'dovich 
1989). The details of the simulation procedure are de- 
scribed in Borgani, Coles & Moscardini (1994) and Pa- 
per I, to which we refer the reader for further details. In 
the following we will sketch only the basic steps of our 
simulation procedure. 

The ZA assumes that, at the time t, the final (Eu- 
lerian) positions x(q,i) of self-gravitating particles are 
given by x(q,i) = q + b(t) V q tp(q). Here q are the ini- 
tial (grid) particle positions, fe(t) is the growing mode of 
the evolution of linear density perturbations and tp(q) is 
the gravitational potential, which is related to the ini- 
tial (linear) density fluctuation field, 6(q), through the 
Poisson equation 'V^tl'(q) = —S(q)/a. Although gravity 
determines the particle displacements according to the 
ZA in a non-local fashion, no gravitational acceleration 
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appears in the equation of motion. Therefore, particles 
move inertially along straight lines. The corresponding 
velocity, 

v = b(t) V„tf( g ) , (13) 

is related to the initial fluctuation field according to lin- 
ear theory. In order to reduce the effects of shell cross- 
ing, we filter out short-wavelength fluctuation modes 
according to P(k) — > P(k) exp( — k 2 R 2 f), where Rf is 
chosen in such a way that the average number of streams 
at each point is 1.1 (see Paper I). 

In this paper we mainly consider simulations of box 
size L = 480 h~ 1 Mpc which are larger than those pre- 
sented in Paper I, which had L = 320 ft. -1 Mpc. Since we 
consider the same number of particles and grid points 
(namely 128 3 ), we have a lower resolution: the grid 
size is 3.75 ft. -1 Mpc. This is not a serious problem for 
the present analysis, since we do not need to know 
the cluster positions with great precision. We merely 
need to check how clusters trace the overall density 
field, and, then, the dynamical origin of the motion in- 
duced by gravity on suitably chosen observers. In fact, 
in Section 3 we have shown that, although the observed 
alignment between the CMB dipole and that traced by 
real clusters already occurs at relatively small distances 
(^ 50 ft. -1 Mpc), the dynamical contribution of clusters 
to the LG motion nevertheless saturates at scales con- 
siderably larger than this (~ 160 ft. -1 Mpc). We therefore 
expect that waves within a box of 320 ft. -1 Mpc cannot 
provide by themselves a completely satisfactory descrip- 
tion of the observational situation. Furthermore, in or- 
der to verify that the adopted box size is large enough, 
we will also refer in some cases to simulations based on 
a larger (L = 640ft _1 Mpc) box size to corroborate our 
results. 

Particles are initially placed on the grid points and 
are moved according to the ZA. The density and velocity 
fields on the grid are then reassigned through a TSC in- 
terpolation scheme (e.g. Hockney & Eastwood 1981) for 
the mass and the momentum carried by each particle. 
Clusters are then identified as the N c i = {Lj {r c i}) 3 high- 
est local density maxima on the grid, where (r c i) is the 
average cluster separation. In the following, we assume 
(Pel) = 38 ft. -1 Mpc, so that we end up with a total of 
2015 clusters in each realization. The above value of (r c i) 
is half-way between that of the Abell and ACO samples 
(~ 40 ft. -1 Mpc and ~ 35 ft. -1 Mpc, respectively). In any 
case, we have verified that no significant changes in the 
statistical properties of the cluster dipole occurs when 
passing from (r c i) = 38 h -1 Mpc to either value of (r c i). 

We ran simulations for the six different models of 
the initial fluctuation spectrum. All the models, except 



the open one, are normalized to be consistent with the 
quadrupole of the CMB temperature anisotropy mea- 
sured by COBE (Bennett et al. 1994). For each model 
we generated 20 random realizations. Four models have 
Q = 1 for the cosmic density parameter: the standard 
CDM model (SCDM), with ft = 0.5 and a s = 1 for the 
r.m.s. fluctuation amplitude within a top-hat sphere of 
8 ft." 1 Mpc; the tilted CDM model (TCDM), with pri- 
mordial spectral index n = 0.7, ft = 0.5 and as = 0.5 
(e.g. Adams et al. 1993; Moscardini et al. 1995); a low 
Hubble constant CDM model (LOWH), with ft = 0.3 
and <j 8 = (1.6)" 1 (Bartlett et al. 1994); the Cold+Hot 
DM model (CHDM), with Q hot = 0.3 for the fractional 
density of the hot component, ft = 0.5 and as = (1.5) -1 
(e.g. Klypin et al. 1993). Other two models have low 
density parameter (f2 = 0.2) both with open (OCDM 
with ft = 1 and as = 1) and with flat geometry (ACDM 
with Qa = 0.8 for the cosmological constant term, ft = 1 
and as = 1.3). A problem with this latter model is that 
the r.m.s. value of as at the cluster scale (~ 8 ft -1 Mpc) 
exceeds unity. We do not therefore expect that the ZA 
will give a correct description of the cluster distribu- 
tion in this case. Although this could be a problem for 
clustering analysis, which requires a precise determina- 
tion of cluster positions, it is of secondary importance 
for the purpose of this work, where we just require that 
clusters roughly trace the large scale DM distribution. 
In a sense, while the determination of cluster correla- 
tions requires the knowledge of cluster positions within 
superclusters, questions concerning the amplitude and 
scale of convergence of the cluster dipole can be prob- 
ably addressed just by knowing the positions of such 
superclusters. 

In a recent paper, Gaztafiaga, Croft & Dalton 
(1995) have attempted to cast doubt upon the relia- 
bility of the use of the ZA for cluster simulations like 
these. Their arguments are, however, based solely on a 
comparison between N-body and ZA approaches for the 
ACDM model, which we already know to be the most 
problematic for our method. In contrast, they fail to dis- 
cuss the fact that, for the SCDM and CHDM models, 
they obtain similar results and arrive at the same con- 
clusions as in our analysis based on the ZA (Plionis et 
al. 1995; Borgani et al. 1995). More details of tests we 
have made demonstrating the general reliability of our 
ZA cluster simulations for clustering studies are given 
in Borgani et al. (1995). 

4.2 LG— observer identification 

In order to perform a consistent comparison between 
simulations and real data, we need to reproduce the ob- 
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servational set-up as closely as possible. To this end, we 
need to define 'observer' characteristics, such as their 
local density and velocity. In particular, a reliable defi- 
nition of the 'observer' peculiar velocity is crucial, for at 
least two reasons: (a) the frequency for the occurrence 
of good alignment between v b s and V c i should tell us 
whether the small value of the observed A9 (cf. Table 
2) is quite typical or it is due to the presence of a highly 
coherent "attractor" in the nearby cluster distribution; 
(b) the comparison of the \v b s \ with \V c i\ could furnish 
an estimate of the parameter (3. 

We place observers at grid points and define both, 
their local density contrast and bulk velocity, averaged 
over a top-hat sphere of radius R = 7.5 ft. -1 Mpc centred 
on them. The local value of S is estimated by convolving 
the Fourier transform of the fluctuation density field 
with the Fourier transform of the top-hat window, 



W(kR) 



3(sin kR — kR cos kR) 
(kR) 3 



(14) 



The definition of the observer's velocity is however less 
straightforward. We have considered two alternative 
ways of defining this quantity: 

• The first is to use eq.(13) to define velocities for the 
particles and then to assign the corresponding mo- 
mentum on the grid by means of the top-hat win- 
dow function. The drawback of this procedure lies 
on the fact that the velocities of eq.(13) are in some 
sense not "synchronized" with the evolution of the 
density field; they refer to the initial field, 6(q), 
while the Zel'dovich evolution has moved structures 
away from their initial locations. Therefore, there is 
no reason to expect a priori any accurate alignment 
between the direction of the linear (initial) velocity 
vector, attached to the final particle position and 
the dipole estimated at that position by integrating 
over the ZA-evolved fluctuation field. 

• The second criterion, which is in fact the one that 
we chose to adopt, is to define the velocity of a given 
'observer' by means of eq.(5), where the dipole term 
is evaluated in a sphere over the whole distribution 
of the 128 3 DM particles at a distance larger than 
R = 7.5 h -1 Mpc from the observer. Note that when 
the borders of the box are encountered we wrap- 
around, so that spherical symmetry is preserved. 
Here the value of the biasing parameter is set to 
unity by definition, while Q is that of the model 
under consideration. This method of assignment as- 
sumes that, after ZA evolution, the gravitational 
potential and velocity fields are still connected by 
the linear relation v cx Vtp, but where tp is now the 



Table 2. Column 2: Number of observers satisfying the 
LG requirements. Column 3: r.m.s. observer velocity (in 
ferns -1 ). Columns 4 and 5: fraction of observers with mis- 
alignment < 20° in LG and CMB frames, respectively. 



Model 


N LG 


^rms 


f{A9 LG ) 


/(A9cmb) 


SCDM 


3019 


592 


0.38 


0.21 


LOWH 


2366 


461 


0.32 


0.20 


TCDM 


1347 


384 


0.28 


0.17 


CHDM 


3192 


593 


0.39 


0.22 


OCDM 


331 


331 


0.41 


0.31 


ACDM 


974 


399 


0.44 


0.32 



potential corresponding to the ZA evolved density 
field. 

Note that we estimate the local S values and the pe- 
culiar velocity for 2000 observers chosen at random po- 
sitions in each realization, so that we end up with a 
total of 4 x 10 4 observers for each DM model. We then 
select those 'observers' that have the same character- 
istics as the Local Group, the number of which can 
also be viewed as a first discriminative test for our DM 
models (cf. Gorsky et al. 1989; Davis, Strauss & Yahil 
1991; Strauss, Cen & Ostriker 1993; Tormen et al. 1993; 
Strauss et al. 1994; Moscardini et al. 1995). These char- 
acteristics are: 

(1) Peculiar velocity Vlg = 627 ± 44km s -1 (error cor- 
responding to 2<j uncertainties; Kogut et al. 1993) for 
a top-hat sphere of radius R = 7.5 h~ 1 Mpc centred 
on the observer; 

(2) Density contrast within the same sphere in the range 
-0.2 < S LG < 1. 

A further requirement concerns the quietness of the lo- 
cal flow, which implies a small value for the local shear. 
Moscardini et al. (1995) argued that this kind of con- 
straint is not very restrictive (see also Tormen et al. 
1993). On the other hand, Schlegel et al. (1994) found 
that CHDM performs much better than SCDM in pro- 
ducing locally quite Hubble flows, which are similar to 
that observed around the Local Group. In any case, we 
do not expect any strong shear to be induced by the 
laminar flow description of the ZA, so we have decided 
not to impose it in our present analysis. 

In Table 2 we report the fraction of observers which 
satisfy the LG requirements for all the models. Also re- 
ported are the r.m.s. values of velocities of these ob- 
servers. In Figure 3 we plot the frequency distributions 
P(S) and P(y ) for the local density and velocities of 
the selected observers. The errorbars in each bin repre- 
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sent the la scatter over the ensemble of 20 realizations, 
which is consequently an estimate of the cosmic vari- 
ance. The solid curves, superimposed on P{y ), are the 
best-fitting Maxwellian distributions. We note that the 
fit is very good in all the cases, showing that the ini- 
tial Gaussian statistics of the velocity field is preserved 
by the Zel'dovich dynamics. The shaded areas delineate 
the observational uncertainties on the S b s and v b s val- 
ues. It is evident that a large fraction of observers for 
all the models has a local density within the observa- 
tional interval. The probability, however, of reproducing 
the LG velocity depends quite strongly on the model. 
For instance, SCDM and CHDM reproduce the obser- 
vational situation with a high probability. On the other 
hand, low-density models give lower velocities, with 
OCDM behaving even worse than ACDM, because of 
the lower normalization of as . Consequently, adopting a 
substantially lower normalization for the OCDM model, 
as — 0.5, according to some proposed COBE normalized 
open models (Kamionkowski & Spergel 1994; Gorski et 
al. 1995), the velocities produced would be even smaller 
and the probability of having a LG-like observer would 
be smaller still. As for TCDM and LOWH models, the 
lower velocities with respect to SCDM are also due to 
their lower as normalization. 

5 RESULTS FROM SIMULATIONS 

In this section we present a dipole analysis of the clus- 
ter simulations. After showing the importance of the 
requirement that CMB and cluster dipoles be aligned, 
we verify the ability of the spherical harmonic recon- 
struction method, which we applied on the real data 
samples, to recover information about the intrinsic (i.e. 
whole-sky) cluster distribution. To this end, we extract 
realistic data sets from the simulations by imposing the 
same observational biases as in the Abell/ACO sample, 
and analyze them on the same footing as in Section 3. 
Finally, we will compare results based on the real data 
and the simulations in order to put constraints on the 
DM models considered. 

5.1 Cluster and CMB dipole alignment 

A fundamental requirement to select realistic observers 
is the existence of a good alignment between the CMB 
and the cluster dipole vectors. Under the hypothesis 
that linear gravitational instability generates the large- 
scale motions and that clusters trace fairly the under- 
lying matter density field, modulo the biasing factor b 
defined in eq.(3), we should expect the two vectors to 



be aligned, as is the case for the observational samples 
(see Table 1). In order to verify whether such a good 
alignment is always to be expected, we estimate, for 
each observer the true cluster dipole vector correspond- 
ing to the whole cluster distribution within 240 ft. -1 Mpc 
(i.e. without including observational selection effects) 
and the misalignment angle A9 with respect to v v . 

In Figure 4 we plot the /3lg values estimated by 
using eq.(5) for 5000 ACDM and CDHM observers as 
a function of the corresponding A9 misalignment an- 
gles between their CMB and cluster dipoles. The dipole 
amplitude is measured for all the observers at the scale 
R-max = 230 ft. -1 Mpc. Also plotted is the correspond- 
ing cumulative A9 distributions, both in the LG (solid 
curves) and in the CMB (dashed curves) frames. At least 
95% of observers have A9 < 90°, thus showing that the 
cluster distribution does indeed reflect the dynamical 
origin of the observers' motion. We can summarize the 
main conclusions of the A9 analysis by noting that: 

• Dipoles in the LG frame have systematically a bet- 
ter alignment than those in the CMB frame, in 
agreement with the qualitative explanation given 
in Section 3.2 and with what we have found in the 
observational case. 

• The difference between LG and CMB frames is 
smaller for ACDM than for CHDM, as expected 
from eq.(12) since the observers have smaller pecu- 
liar velocities in the former model. 

• Due probably to the sparseness of the cluster distri- 
bution, most simulation observers measure a much 
larger value of A9 than the values found for the 
Abell/ACO clusters and there is therefore no a pri- 
ori reason to expect an alignment as accurate as 
that seen in the real data. This shows that the LG 
motion is probably generated by a large, very coher- 
ent structure, which is well sampled by the observed 
cluster distribution. 

In order to reproduce the observational situation in the 
following analysis, we will select as LG candidates only 
those observers that measure A9 < 20° in both CMB 
and LG frames. This is quite a conservative restriction 
since, as shown in Figure 4, no substantial changes in 
the (3 estimates should occur for A0< 40°. In fact, we 
have verified that taking this larger limit for A9 does 
not change results in either the LG or CMB frame. In 
Table 2 we list the percentage of observers passing only 
this test for the different models. 

Although no stringent constraints can be posed 
on the models from the statistics of A9, it turns out 
that low-density models produce more relatively bet- 
ter "aligned" observers, especially in the CMB frame. 
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This could well be due to the fact that such models, 
which have more large-scale power, produce more co- 
herent and extended structures in the cluster distribu- 
tion. In the cases with Q = 1 observers, the less accu- 
rate alignment in the CMB frame is compensated in the 
LG frame by their higher velocities with respect to the 
Q = 0.2 models. 



5.2 Estimates of the (3— parameter 

5.2.1 Whole box analysis 

In Figure 5 we plot the frequency distribution of (3 val- 
ues measured by all the observers passing the A9 test, 
for both ACDM and CHDM. In each panel, we report 
results of the analysis done in real space (solid line), 
in the CMB frame (dotted line) and in the LG frame 
(dashed line). We conclude that: 

• The (3 estimates in the CMB frame are always very 
close to the "true" ones of real space, which is ev- 
idence that most LG-like observers in our simula- 
tions are not embedded in any coherent bulk flow 
(see discussion in Section 3.2)*. 

• The (3 estimates in the LG frame are systematically 
larger than those in the CMB frame, which is to be 
expected according to eq.(12). 

• The ACDM model produces systematically lower 
(3 values than CHDM. The P(f3) distributions are, 
however, too broad to allow strong discrimination 
even between these two cases. 

In order to verify whether the presence of fluctuation 
modes on scales larger than the box size can alter our 
results, we also analyzed simulations in a box of side 
L = 640 ft. -1 Mpc. We verify that increasing the box size 
does not alter the P(f3) distributions significantly, while 
we also obtain very similar values of 0. For example, 
in the ACDM model, we get j3 C MB = 0.16 ± 0.06 for 
the L = 480/i _1 Mpc box and ficMB = 0.16 zb 0.08 for 
the L = 640 h -1 Mpc box. Therefore, wavelength modes 
on scales comparable to or exceeding the box size used, 
act in the same way on the value of V c i and on the 
observer velocity v p , so as to make their ratio (i.e. the 
f3 parameter) almost independent of the box size. 



* In future we plan to investigate whether the existence or 
not of 'coherent motions' in our simulations could discrimi- 
nate between the different DM models. 



Table 3. Columns 2 and 4: Average /3 value for the mock 
samples, after correction for the 0(#) mask in the CMB and 
LG frames, respectively. Columns 3 and 5: Number of LG- 
like observers measuring a (3 value in the observational range 
after correcting for the 0(#) mask for CMB and LG analyses, 
respectively. 



Model CMB frame LG frame 







Mask 


iV real 


Mask 


iV real 


SCDM 


0. 


,23 ± 0.14 


181 


0.15 ±0.11 


492 


LOWH 


0. 


,20 ±0.12 


107 


0.13 ±0.15 


195 


TCDM 


0. 


,19 ± 0.11 


170 


0.12 ±0.09 


340 


CHDM 


0. 


,22 ± 0.13 


238 


0.14 ±0.10 


435 


OCDM 


0. 


,16 ± 0.07 


45 


0.09 ±0.06 


75 


ACDM 


0. 


,17 ±0.09 


122 


0.11 ±0.07 


216 



5.2.2 Analysis after including realistic selection 
effects 

Observational biases, like Galactic obscuration and red- 
shift selection, are expected to affect dipole determi- 
nations, especially when they are not isotropically dis- 
tributed around the sky. Furthermore, the redshift se- 
lection function could significantly affect the dipole es- 
timation especially if the redshift completeness limit is 
lower than the true dipole convergence depth. However, 
since we limit our analysis to R ma x = 230 ft. -1 Mpc, a 
depth at which the cluster dipoles have converged to 
their final values and at which the observed Abell/ACO 
cluster sample is quite complete, its effect should be 
negligible. 

In order to account properly for these effects, we 
extract cluster samples with the same 0($) and P(z) 
selection functions as the Abell/ACO sample [see eqs. 
(11) and (10), respectively], for each LG-like observer 
passing the misalignment test. Since the effect of the 
0($) mask could depend on the Galactic latitude of 
the CMB dipole direction, we choose the observer co- 
ordinate system so that its v v has the same direction 
(I = 277°, b = 30°) as the observed CMB dipole. For 
each of these observers we construct mock catalogues 
by performing 100 Monte-Carlo realizations of the se- 
lection functions, in order to minimize the noise due to 
accidental cluster selection in regions of strong Galactic 
obscuration. The resulting dipole amplitude of a given 
observer is taken to be the average over all those Monte 
Carlo realizations for which the misalignment angle still 
satisfies the constraint A9 < 20° after the mask correc- 
tion has been applied (as in the real cluster case). 

In Figure 6 we plot the frequency distribution of the 
(3 values estimated in the CMB frame, for all the LG-like 
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ACDM and CHDM observers, passing the three tests for 
vlg, &lg and A9. Intrinsic dipole /3-values (solid his- 
tograms) are compared to those before mask correction 
(left panels) and after mask correction (right panels). 
The effect of introducing observational biases is that of 
increasing the noise, i.e. the width of the distribution. 
As far as the mask correction procedure is concerned we 
can conclude from this study that: 

• Dipole alignment: The number of well-aligned 
dipoles after correcting for Galactic absorption is 
~ 20 — 30% larger than that before the corrections; 
this result shows that the correction procedure is 
effective in recovering on average the cluster dipole 
direction. 

• Dipole amplitude: The intrinsic P(f3) distribution 
is recovered quite well after the mask correction, 
confirming again the reliability of our procedure. 

In Figure 6 we also plot, as the dashed band, the la 
uncertainty in the the Abell/ACO (3 estimate (c.f. Ta- 
ble 1). We note that the observational /3-range is more 
typical in a CHDM universe, while it falls on the high-/3 
tail of the P((3) distribution for the ACDM model. 

In Table 3 we report the average values, /3, mea- 
sured by LG-like observers for the realistic cluster sam- 
ples, in both the LG and CMB frames, for all the models 
(after correcting for Galactic absorbtion). Although not 
reported here, the whole-box values are very similar 
to those of Table 3, as indicated by the similarity be- 
tween the corresponding P(f3) distributions (see Figure 
6). Also reported is the number N rea i of such observers, 
which besides satisfying the LG and alignment tests, 
measure a (3 value consistent with the real one. 

Despite the fact that differences exist between the 
different models, all of them produce an appreciable 
number of such observers. This means that: 

• It is quite difficult to discriminate at a high confi- 
dence level between different initial power-spectra, 
which is supported also by the large uncertainties 
in the (3 estimate due to the large spread of the 
P(f3) distributions. 

• The Abell/ACO cluster dipole cannot be taken as 
strong evidence for a low-f2 Universe, as has been 
suggested by Scaramella et al. (1991) and PV91. 
Instead, the indications are, if anything, that some 
n o = 1 models are marginally preferred to low- 
density (f2 = 0.2) models. 

It is interesting to note that the Abell/ACO data show a 
difference between /3lg and /3cmb which is larger than 
the typical difference found in the simulated data sets. 
This further indicates that the observational situation is 



probably characterized by the presence of a large coher- 
ent structure, which causes a bulk motion of the Local 
Universe of an amplitude similar to the velocity of the 
Local Group (cf. PV91; Plionis, Coles & Catelan 1993; 
Branchini & Plionis 1995) thus generating, in the man- 
ner described by eq.(12), relatively larger differences be- 
tween LG- and CMB-estimated cluster distances than 
those found in the simulations. 



5.3 Estimating the cosmological density 
parameter 

We can now investigate whether the estimated values 
of (3 turn into a reliable estimate of the density pa- 
rameter Q . In order to pass from (3 to Q one must 
know, in principle, the local biasing parameter, i.e. the 
proportionality constant which relates cluster and DM 
distributions point by point. What is usually estimated, 
however, is the global biasing factor b c i, evaluated in 
a statistical sense from global correlation properties of 
the cluster distribution. Nonetheless, the non-linear se- 
lection of clusters as peaks of the density field means 
that the correct biasing parameter is actually a func- 
tion of the position. 

In the simulations we estimate b c i statistically, as 
the ratio between the r.m.s. fluctuation amplitude for 
cluster and DM distributions at different scales, 



b c i{R) 



o-d(R) 



(15) 



o~dm{R) 

We estimate the r.m.s. fluctuations in the DM distribu- 
tion as 



°~dm{R) 



— f 



dk k 2 P(k)W 2 (kR), 



(16) 



where W{kK) is the Fourier transform of the top-hat 
window function, as given by eq.(14), and P(k) the lin- 
ear power-spectrum. 

Since the dipole analysis is based on counts within 
spheres of different radii, we estimate <r c i(i?) through 
counts of clusters within top-hat spheres of radius R. 
Furthermore, as in the cluster dipole analysis, we do 
not correct the variance of the cluster distribution for 
Poissonian shot-noise. The analysis is done by ran- 
domly placing 5000 spheres within each simulation box, 
with radii ranging from 50 h~ 1 Mpc and 120 ft. -1 Mpc. 
Although we need to estimate the value of b c i at the 
largest considered scale, R ma x = 230 ft. -1 Mpc, it turns 
out that, at this scale, the variance analysis is domi- 
nated by the finite box size by the consequent effect 
of periodic boundary conditions. However, consistently 
with what we found in Paper I, we have verified that 
b c i is almost independent of scale in the range we have 
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Table 4. Column 2: Cluster biasing parameters and relative 
"cosmic" r.m.s. uncertainties. Column 3: fi estimates for 
LG— like observers in the CMB frame for whole— box cluster 
samples. Column 4: The same as in column 3, but after in- 
troducing the 0(#) mask and correcting for it. Errors come 
from uncertainties in the (3 estimate (cf. Table 3). 



Model 



(f} }(whole sky) {f} )(masked) 



SCDM 


3.74 ±0.17 


1 o +1 - 3 

I - u -0.8 


o.8; 


LOWH 


5.09 ±0.22 


1 1 + 12 


l.o; 


TCDM 


7.43 ±0.39 


1 9+ 21 


1.8; 


CHDM 


4.11 ±0.16 


1 0+ 13 


0.9; 


OCDM 


2.96 ±0.14 


o.23l° ;J 5 3 


0.29 


ACDM 


2.36 ±0.11 


0.20t° ;» 


0.22 



+1.1 

-0.6 
+ 1.3 
-0.8 
+ 2.0 
-1.4 
+ 1.5 
-0.7 
+ 0.24 
-0.18 
+ 0.22 
-0.16 



explored. We therefore feel justified in extrapolating 
the value of b c i determined on smaller scales up to 
R-max = 230 ft. -1 Mpc. The b c i values for all six of our 
models are reported in Table 4, along with the corre- 
sponding "cosmic" r.m.s. uncertainties, estimated from 
the 20 realizations of each model. The central value is 
used to convert the (3 estimate of each observer into the 
corresponding Q estimate. 

In Figure 7 we plot the distribution of Q estimated 
in the CMB frame for the CHDM and ACDM LG- 
like observers. For both models, the P(Q ) distribution 
turns out to be rather broad and a non-negligible prob- 
ability exists to measure values of Q substantially dif- 
ferent than the correct ones. In Table 4 we report, for all 
the models, the average Q values measured by LG-like 
observers in the CMB frame. The two columns refer to 
estimates from the whole-box cluster distributions and 
from the mock cluster samples. Although from Figure 7 
we note that the P(Q ) peaks at Q values smaller than 
the correct ones, the average value always agrees well 
with the "intrinsic" Q . As for the comparison between 
low-density and high-density models, it turns out that 
there is a non-negligible probability of measuring Q as 
low as 0.2 in a Q = 1 universe. Vice versa, if Q = 0.2 it 
is almost impossible to measure, by chance, Q = 1. On 
the other hand, values as large as 0.4 are commonplace 
is both of these models. 



6 DISCUSSION AND CONCLUSIONS 

We have compared the 3D dipole structure of an ex- 
tended sample of Abell/ACO clusters, supplemented 
with the new redshifts of Katgert et al. (1995), to that 
generated by numerical simulations of the cluster distri- 
bution, based on six different DM models. The simula- 



tion technique, based on the Zel'dovich approximation, 
is so computationally cheap that it allows us to run a 
set of 20 realizations for each model. We present our 
main conclusions in the following subsections: 

6.1 Real Cluster Dipole Results 

We have estimated the Abell/ACO cluster dipole ac- 
counting for the effects of Galactic obscuration by ap- 
plying the spherical harmonic method (see Appendix). 
Our results (see Table 1 and Figure 2) are in general 
agreement with those coming from other studies, based 
however on fewer measured cluster redshifts (PV91; 
Scaramella et al. 1991; Branchini & Plionis 1995) and 
can be summarized as follows. 

(a) The dipole convergence occurs at a scale, R CO nv — 
160 ft. -1 Mpc, substantially larger than the conver- 
gence scale indicated by shallower galaxy samples 
(e.g. Strauss & Davis 1988; Lynden-Bell et al. 1989; 
Rowan-Robinson et al. 1990), but still well within 
the completeness limit (~ 240 h -1 Mpc) of the cluster 
sample. 

(b) Despite the large value of R CO nv , the direction of the 
cluster dipole is already very well aligned (A0< 20°) 
with the CMB dipole (Kogut et al. 1993) at much 
smaller scales (~ 50 ft. -1 Mpc). This suggests the pres- 
ence of an extended coherent puller generating the 
motion of the Local Group. 

(c) The final value of the dipole amplitude, estimated 
both in LG and CMB frames, turns into values of 
the /3-parameter of /3lg = 0.15 ± 0.04 and /3cmb = 
0.25 ± 0.06. The difference between these two values 
is due to the different distances at which clusters are 
placed, according to eq.(12), along the direction of 
our motion. 

6.2 Simulation Cluster Dipole Results 

We took particular care to reproduce the actual obser- 
vational setup as closely as possible in our simulated 
cluster distributions. This was realized in two steps: 

• We identify those 'observers' which have local (i.e. 
within a sphere of 7.5 h -1 Mpc radius) values for the 
density fluctuation S b s and peculiar velocity v b s 
consistent with those of the LG; i.e. —0.2 < Slg < 1 
and vlg = 627 ± 44 km s -1 . As a further criterion, 
we also require that the misalignment angle A9 be- 
tween the direction of the whole-sky cluster dipole 
and the observer velocity is A9 < 20°. 

• For each of these LG-like observers we generate 
mock cluster samples having the same character- 
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istics (i.e. number of clusters, Galactic latitude and 
redshift selection functions) as the real sample. 
Also, for each observer the coordinate system is 
suitably rotated so that the vector \ bs points in 
the same direction as the observed CMB dipole. 

The main results concerning the simulation analysis 
and their comparison with real data can be summarized 
as follows. 

(a) The LG-like observers found in each model (Table 
2 and Figure 3) are quite typical in Q = 1 models 
with a moderate as normalization, like SCDM and 
CHDM, but are much less likely in the Q = 0.2 mod- 
els (OCDM and ACDM) and in a fi = 1 models with 
a low normalization [as — 0.5), like TCDM. 

(b) Regarding the frequency of occurrence of small mis- 
alignment angles between the cluster dipoles and the 
observer velocities, we find that small A9 values are 
not to be expected in most cases. In fact, A9 < 20° 
is measured by 30-40% of the observers in LG frame 
and by 20-30% of them in CMB frame (Table 2 and 
Figure 4). 

(c) The (3 estimates for real-space and CMB frame are 
quite similar, while analyses in the LG frame give sys- 
tematically lower /3's. Although we estimate /3lg < 
/3c MB also for real Abell/ACO clusters, Branchini & 
Plionis (1995) found that their real-space estimate of 
(3 is about half way between /3lg and /3cmb- This 
indicates that the observational situation is charac- 
terized by a coherent motion of the Local Group, 
which also involves the nearest clusters (see also Ta- 
ble 7 of PV91 and Plionis, Coles & Catelan 1993). In 
fact, if this is the case, then eq.(12) would imply that 
the cluster redshifts in the CMB frame are artificially 
overestimated. 

(d) Although the P(f3) distributions are significantly dif- 
ferent for different models, especially when comparing 
low-density and high-density models (see Figure 6), 
both the width of P(f3) and the rather large uncer- 
tainties of the observational (3 value make it difficult 
to constrain any model. We point out, however, that 
SCDM and CHDM models are marginally preferred 
over the others, in the sense that they have the largest 
number of observers measuring a (3 value consistent 
with the real one, while LOWH, OCDM and ACDM 
are less favoured (see Table 3). 

(e) By directly estimating the cluster biasing parameter 
b c i in our simulations, we can determine the value of 
Q that would be estimated by each observer. It turns 
out that single observers can measure Q values very 
different from the true value. The correct value is, 



however, nicely recovered on average (see Table 4 and 
Figure 7). 

6.3 Final Comments 

Based on these results, we can safely state that the an- 
swers to the two questions raised in the Introduction are 
both negative. That is: 

• Due (probably) to the sparse sampling of the under- 
lying density field and the effect of cosmic variance, 
the estimated value of (3 from Abell/ ACO type clus- 
ter distributions can be very different from the true 
value. It therefore seems that we cannot infer the 
value of Q with precision. 

• No strong constraints can be posed on DM models 
based solely on the study of the dipole structure 
(its large amplitude and convergence depth) of the 
cluster distribution. 

A further possibility to discriminate between low- and 
high-density cosmological models is offered by study of 
the large-scale cluster peculiar velocity field (see, e.g., 
Bahcall, Gramann & Cen 1994; Cen, Bahcall & Gra- 
mann 1994; Croft & Efstathiou 1994). The recent re- 
sults of Lauer & Postman (1994) on the presence of 
a large amplitude cluster bulk flow, within ^ 15000 
km/sec, presents a strong challenge for current cosmo- 
logical models (Strauss et al. 1994; Feldman & Watkins 
1994; Jaffe & Kaiser 1994). However, Branchini & Plio- 
nis (1995) and Branchini, Plionis & Sciama (1995), us- 
ing a reconstruction method and linear gravitational 
instability theory to obtain the 3D cluster positions 
and peculiar velocities, found that the Lauer & Post- 
man cluster sample exhibits a bulk flow of only ~ 200 
km/sec, in much better agreement with most cosmolog- 
ical models. 

We plan in the near future to take advantage of 
our ZA-based cluster simulations to test several DM 
spectra against observations of large-scale cluster mo- 
tions, with particular emphasis on distinguishing be- 
tween high- and low-density models. 
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APPENDIX Al: SPHERICAL HARMONICS 
MASK METHOD 

The method we use to treat the Galactic absorption 
effects is based on expanding the sky surface density 
field <r($, ip) in spherical harmonics according to 



(Al) 



1 = 

I 



+ ^ ^ P; m (cos $) (af 1 cos mip + fej 11 sin mip) , 

m= 1 

where < # < 7r, < ip < 2-k and _P, m (cosi?) are 
the Legendre polinomials. Accordingly, the monopole is 
defined as M = a° , while |D| = [(al) 2 + {a{) 2 + {b{) 2 ] 1/2 
for the dipole amplitude. 

If E($, ip) is the observed surface density field and 
<r($, ip) is the intrinsic one, the two are related by the 
mask function Q(D,tp), which defines the missing parts 
of the survey, according to 



(A2) 



If we are interested in recovering the dipole (I = 1) com- 
ponents of E($, ip), correction terms in 0($, <p) should at 
least involve the quadrupole (I = 2) terms. We show be- 
low the detailed procedure to calculate the (a™ , b™) co- 
efficients for the spherical harmonic expansion of <r($, ip) 
in terms of the observed (A™,!)™) coefficients, for a 
given mask model Q(D,tp). Since in the case of Galac- 
tic absorption there is no ip dependence of the mask we 
have that Q(#,<p) = 0(#). 

Starting from the expansion of E($, ip) up to the 
quadrupole order and allowing for the orthogonality re- 
lation of the Legendre polynomials, 

2S W (l + m)\ 



L 



pr{p)p?{p)d{p) 



(A3) 



(2Z + l)(Z-m)! 

(m < 1,1'), we can express the observed coefficients in 
terms of the intrinsic ones as 
Am _ (2l + l)(l-m)\ 



2tt 

»2tt 



(A4) 



(l + m)\ 

I I f^)0(i9)P i m (cos cos rrwp sin fidfidip 

Jo Jo 



where % = 7r/2 — |&|i»m and we introduce 
Hi = H n {A N ) + H n (A s ) 
H n = H^An) — H n (As) , 



(AT) 



where An and As are the values of the Galactic extinc- 
tion amplitude in the North and South Galactic Caps, 
respectively. Therefore, substituting the expansion (Al) 
of <r($, ip) and the above expression for the mask (11), 
eqs.(A4) and (A5) give a linear system of 9 equations 
with the 9 unknown quantities a™, ft™ (' ^ 2, m < I) 
which can be written as: 



a, 



o 



o 

3a, 



(ZH~ - H~) 



4 
4 

5dn 



361 (H+-Ht) + *±(Hr-H 6 -) 



4 
4 

5a? 



{ZHi-Hi) + ^{ZH--H-) + 



5a, 



Al 



Bl 



( Hr -H-) + 1 -^ {Hi -Hi) 
5bl ( Hr -H-) + 1 -^ (Hi -Hi) 



4 

5&J 
4 

15o| 

16 
156j 

16 



(H+-2Hi+Hi) 
(H+ -2Hi +Hi) . 



By inverting the above system of equations we can work 
out the corrected monopole, dipole and quadrupole co- 
efficients in terms of the observed ones. 



BT 



(2Z + 1) (Z - m)! 
2-k (l + m)\ 



(A5) 



y»27T y»7T 

/ / <t($, (p)0($)P; m (cos $) sin mip sin $d$d<p . 
Jo Jo 

We can express the mask 0($) as a function of the dif- 
ferent Galactic absorption coefficients (for the NGC and 
SGC) by defining the following quantity: 

-A 



H n (A) 



J cos 
Jo 



-& sin -& exp 



cos D\ 



dd , (A6) 
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Figure Captions 

Figure 1. The Abell (solid line) and ACO (dashed line) 
cluster density for \b\ > 30°, as a function of redshift. 
Poissonian errorbars are shown only for the Abell sam- 
ple. 

Figure 2. The dipole amplitude V c i for the Abell/ACO 
cluster sample, estimated both in the CMB (open cir- 
cles) and in the LG (filled circles) frames after the spher- 
ical harmonic correction for the effects of Galactic ex- 
tinction. The errorbars at the largest considered scale 
include both Poissonian errors and the effects of the un- 
certainties on the shape of the selection functions (see 
text). 

Figure 3. The frequency distribution for the local den- 
sity contrast S b s and the peculiar velocity v b s within 
a top-hat sphere of 7.5 h~ 1 Mpc radius, for the 40,000 
observers selected in each model. Errorbars correspond 
to the scatter between the 20 realizations of each model, 
while the solid curves for P{y ) are the Maxwellian best- 
fit. Shaded areas correspond to the observational val- 
ues for local density S b s and velocity v b s of the Local 
Group (see text). 

Figure 4. Scatter plots of the /3-parameters, measured 
in LG frame by 5000 observers in ACDM (left panel) 
and CHDM (right panel) simulations, as a function of 
the misalignment angle A9 between the directions of 
CMB and cluster dipole. Results refer to the whole-box 
simulated cluster distributions seen by each observer. 
Also shown in the top panels is the probability of having 
an alignment better that a given A9, both in LG (solid 
curves) and in CMB (dashed curves) frames. 
Figure 5. The /3-parameter frequency distributions as 
measured by LG-like observers, identified in ACDM 
(left panel) and CHDM (right panel) simulations. The 
results come from the whole-box cluster distributions. 
Solid, dotted and dashed histograms refer to analyses 
realized in real space, CMB frame and LG frame, re- 
spectively. 

Figure 6. The frequency distributions of (3 values mea- 
sured in the CMB frame by LG-like observers identified 
in CHDM (upper panels) and in ACDM (lower panels) 
simulations. Solid histograms refer to the whole-box 
cluster distributions while dashed histograms refer to 
the analysis of the mock cluster samples. Left and right 
panels refer to the sample analysis before and after cor- 
recting for Galactic absorbtion mask, respectively. The 
shaded areas correspond to the la interval in the esti- 
mate of the (3 parameter for the real Abell/ACO sample. 
Figure 7. The frequency distributions of the Q val- 
ues measured by LG-like observers in the ACDM (left 
panel) and CHDM (right panel) simulations. 
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